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Forces on Cylinders and Plates in an Oscillating Fluid ! 

Garbis H. Keulegan and Lloyd H. Carpenter 

The inertia and drag coefficients of cylinders and plates in simple sinusoidal currents 
are investigated. The midsection of a rectangular basin with standing waves surging in it is 
selected as the locale of currents. The cylinders and plates are fixed horizontally and below 
the water surface. The average values of the inertia and drag coefficients over a wave cycle 
show variations when the intensity of the current and the size of the cylinders or plates are 
changed. These variations, however, can be correlated with the period parameter U m T/D, 
where U m is the maximum intensity of the sinusoidal current, T is the period of the wave I 

and D is the diameter of the cylinder or the width of the plate. For the cylinders U m T/D 
equaling 15 is a critical condition yielding the lowest value of the inertia coefficient and the 
largest value of the drag coefficient. Tor the plates the higher values of the drag coefficient 
are associated with the smaller values of U m T/D and the higher values of the mass coefficient 
with the larger values of U m T/D. The variation of the coefficients with the phase of the 
wave is examined and the bearing of this on the formula for the forces is discussed. The flow 
patterns around the cylinders and plates are examined photographically, and a suggestion is 
advanced as to the physical meaning of the parameter U m T/D. 



1. Introduction 

In a remarkable paper on the motion of pendulums 
Stokes showed that the expression for the force on a 
sphere oscillating- in an unlimited viscous fluid con- 
sists of two terms, one involving the acceleration of 
the sphere and the other the velocity [l]. 2 Further- 
more, the inertia coefficient involved in the accelera- 
tion term is modified because of viscosity and, 
indeed, is augmented over the theoretical value valid 
for irrotational flow. The drag coefficient associated 
with the velocity term is modified because of the 
acceleration, and its value is greater than it would 
be if the sphere were moving with a constant velocity. 
Subsequent to Stokes' studies, the forces on a sphere 
moving in a viscous fluid in an arbitrary manner 
were investigated by Boussinesq and also by Basset 
[2, 3]. They found that the force experienced by a 
sphere at a given time depends, in general, on the 
entire history of its acceleration as well as the instan- 
taneous velocity and acceleration. As an example, 
if a sphere is accelerated, say with a constant accel- 
eration, from a position of rest to a finite velocity 
and is then kept at this velocity, the force during 
the initial instants of uniform velocity differs from 
the force occurring at a later time. Rayleigh has 
given the formula for the force for this case [4]. The 
force expression of Boussinesq-Basset contains three 
terms, one of which is in the form of an integral 
involving the history of acceleration. If the integral 
is evaluated when the acceleration is represented by 
a sinusoidal function it then yields the modifications 
of the inertia and drag coefficients in Stokes' formula. 

One expects quantitatively different results when 
the oscillating velocities are large and the flow 
turbulent. As yet a theoretical analysis of the 
problem is difficult and much of the desired informa- 
tion must be obtained experimentally. In this 
respect the experimental studies have been dealt 
with variously. One method is due to McNown 
and Wolf [5], who considered the force on a two- 
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dimensional object immersed 
of three parts: 



flow as made up 



F=A 



oP 



d(kU) 
dt 



-fp x dS+~C d D P U\Ul (1) 



where F is the force per unit length in the direction 
of flow, X) U the velocity at points far removed from 
the object; p x , the ^-component of the ambient 
pressure in the absence of tin 1 body; dS, an element 
of the surface area; (\,, the coefficient of drag; and 
k, the virtual mass coefficient. The dimension of the 
body normal to the How is I), and A {) is a circular 
area, A =tD 2 /4:, to which the added mass is referred. 
If A is the cross-sectional area of the body, A = rA , 
r being a ratio, then 



f- 



p z dS=prA 



and finally 



F-A+[m 



dU 

It' 



1 



dt 



■r^]+^C*D P U\U\. 



In this approach the variability of the mass co- 
efficient, k, is implied. Thus, introducing a new 
coefficient k' such that 

and putting 

C m =(k'+r), (3) 

there is obtained from eq (1), the expression 

F=C mP A ^+i CaDpUlUl, (4) 

which in fact constitutes a second approach utilized 
first by Morison and coinvestigators [6, 7]. The 
form of the expression is in agreement with the 
Stokes formula for force on a sphere oscillating in a 
viscous medium. In a general sense one may still 
regard C m as a kind of mass or inertia coefficient. 
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A third approach was proposed by Iversen and 
Balent, who considered the force on an accelerated 
disk moving in one direction [8]. Briefly, 



where 



F=C P DU 2 , 
n _ n (DU D dU\ 

L - L \^rw~di)' 



(5) 



Keim has considered the case of accelerated cylinders 
[9] and Bugliarello that of accelerated spheres [10], 
all motions being in one direction. Here the resort 
it to a single coefficient Q and attempts to separate 
the effects of acceleration and viscosity have not 
been shown to be successful. Accordingly, the 
adoption of this method can have a meaning only 
for mono tonic motions subject to definite limitations 
as to initial and final conditions. 

For oscillatory motions, although the forces are 
more accurately described either using eq (2) or 
eq (4), the latter might be preferred provided the 
coefficients C m and C d could be predicted with some 
precision. The application of the expression to 
vertical piling and large submerged objects by Reid 
and Bretschneider stresses the necessity of having 
these coefficients better determined [11]. 

On the basis of irrotational flow around the cyl- 
inder, C m should equal 2, and one may suppose that 
the value of C d should be identical with that appli- 
cable to a constant velocity. Morison and coinvesti- 
gators have obtained the values of C d and C m in 
particular cases by considering the observed forces 
in the phases of the wave cycle where U or dU/dt 
vanishes. Such determinations show considerable 
variations of C m from the theoretical value and of 
C d from the steady state value at the corresponding 
Reynolds number. Dealing with field studies at 
Caplen, Texas, R. O. Reid found similar variations 
in C m and C d [12]. The variations in the coefficients, 
however, have not yet been correlated with any 
appropriate parameter. 

The present investigation was undertaken with the 
following two objectives in mind. The first was in 
regard to a supplementary function AR that could 
be introduced in eq (4) for a truer representation of 
force when considering the coefficients C m and C d as 
being constant throughout a given wave cycle. The 
necessity for the term AR is associated with the 
eventuality that the point values of C m and C d 
deviate from their average values. The second 
objective was to examine the possibility of correlating 
the average values of C m and C d with a parameter 
U m T/D, where U m is the amplitude of the harmoni- 
cally varying velocity, T is the period of the oscilla- 
tions, and D is the diameter of a cylinder or the 
breadth of a rectangular plate. The mid-cross 
section of a large rectangular vessel with standing 
waves surging in it was chosen as the field of harmoni- 
cally varying current. The cylinders and plates were 
held fixed horizontally, totally submerged in water 
and extending from one side of the vessel to the other 
to approximate as closely as possible the condition 
of infinite length. 



2. Fluid Forces on an Immersed Body at 
Rest in a Moving Liquid 

It would be instructive to consider the momentum 
equations discussed by Murnaghan for the evaluation 
of force on objects immersed in a perfect liquid [13]. 
The method, however, is now generalized to apply 
to imperfect liquids. 

Consider the case of two-dimensional flow with x 
horizontal and z vertical. The equation of motion 
in the ^-direction is 



p \dt +U dx +W dz)-^x r+ ^ 



(6) 



where u and w are the velocity components along 
the axes x and z, p the density of the liquid, p xx the 
normal stress on an elementary surface perpendicular 
to x, and p zx the tangential stress on an elementary 
surface normal to z, the stress being in the direction 
of x. Because of the incompressibility of the liquid, 



bu dw 
dx dz 



and eq (2) becomes 



du / d 2 . d 



bz ) dx 



dx 



(7) 



(8) 



Take the immersed cylindrical body of surface S, 
as in figure 1, and draw a surface S' of arbitrary shape 
which encloses the cylinder. Let w be the region 
bounded by S and S' and I and n the direction cosines 
of the normal drawn inward into the region. Inte- 
grating eq (8) throughout co, and in this making use of 
Green's Theorem, one finds 

p ^-jdu—p u(lu J rnw)dS—p u(lu+nw)dS'= 

—J (lPxz+np zx )dS— I (lp xz +np 2X )dS\ (9) 





Figure 1. Notation diagram for force analysis. 
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Over the surface S of the immersed body lu-\-nw van- 
ishes because the body is at rest. Also I (lp xx -{- 

np zx )dS=F, that is, the ^-component of the force ex- 
erted on the solid by the moving liquid. It may be 
assumed that if S' is removed sufficiently from the 
body the tangential stress p zx on S' vanishes and the 
normal stress p xx reduces to the hydrostatic pressure 
— p. Solving for F, 



F= 



-p f |f^>+P i u(lu+nw)dS'+ flpdS'. (10) 



The later relation may be given in another form, 
suitable for the present purpose. Select the bound- 
ing surface S' as the rectangular strip shown in 
figure 1. The plane S[ to the left of the cylinder 
passes through the point %— —X\ and the plane S 2 to 
the right passes through x=x x . Denoting the hori- 
zontal velocity components at the points Pi and P 2 
with the comxnon elevation z x by u x and u 2 , and the 
pressures by />, and p 2 , eq (10) now reduces to 

F= —p I jfid<*+p I (u\— u 2 2 )dz 1 + I (pi—p^dzi, 



(ID 



which is the momentum equation of familiar form. 

This may be specialized to evaluate the force on a 
circular cylinder when the motion is irrotational. 
Letting U be the undisturbed velocity and referring 
to Lamb [14], 



*=d~l+^cos20l 



w=-U-, sin 20 
■pi 



p dU ( . a 2 \ . . 1 , « . ^ 
- = -77 ( H — ) sin 0— - (u 2 +w 2 ) 
p at \ r / 2 



(12) 



where a is the radius of the cylinder, r is radial dis- 
tance, and is the angle between a radius vector and 
the vertical line x=0 passing through the center of 
the cylinder. Clearly, Ui=u 2 and the momentum 
equation, eq (11) , reduces to 

Introducing the values of u and p from eq (12), and 
omitting the straightforward but somewhat lengthy 
evaluations, the result is 



*=2,fa.p, 



or in terms of the diameter D of the cylinder 



Fx=C n 



pwD 2 dU 
4 dt' 



(14) 



where C m =2. 

Next, suppose that the undisturbed velocity is con- 
stant and that the body experiences a drag. With the 
liquid extending to infinity and ignoring the variation 
of pressures from the shedding eddies, or, more prop- 
erly, assuming that the surfaces 8[ and #2 are far 
removed from the cylinder, p x =p 2l and eq (11) re- 
duces to 



F 2 = 



■r. 



(vl—v%)dzi. 



(15) 



The velocity Ui=U, and u 2 =mU, where m is de- 
pendent on zJD and on Reynolds number UD/v. 
Thus, 

F 2 =C d pD^- (16) 



where 



a=2 






(1 — ?n 2 )d 



D 



It appears that in ordinary cases where the flow 
departs from irrotat ionalitv and becomes unsteady 
and eddying, eq (11) is still the basis for evaluating 
the force, since the first and third integrals may be 
associated with acceleration and the second with 
drag. That is, the coefficients C m and C d are de- 
rived from eq (13) and (15) pro vie led the velocities 
and pressures can be given. The force of the state- 
ment is only academic, since in the (lows involving 
separation and intermittent eddy formation the 
pressures and velocities are not known and the 
integrations in eq (11), at present, cannot be carried 
out. Nevertheless, experience suggests that eq 
(4) remains useful at least for sinusoidal motions, if 
allowance can be made for the variations in C m 
and Cd. 

Had one carried out the integrations in eq (11) 
for an extended plate using the known velocity ex- 
pressions derived from the Kirchoff solution for the 
impact on a lamina, definite values for C m and C d 
would have resulted. This would have shown in 
principle the existence of a relation between C m 
and C d in the absence of eddy formation. In the 
Kirchoff solution the wake is of infinite length and 
this is cause for concern. McNown overcomes this 
difficulty by considering the case of a closed wake 
as between two plates and finds a relation between 
k and C d or between C m and C d [15]. This result is 
very significant as it points to the path to be fol- 
lowed in analytical approaches taking into account 
also the effect of the eddy processes. With cylinders 
the changing separation seats are a cause of added 
difficulty. 

Meanwhile, the tasks of the experimental investi- 
gations become more necessary. Not only are the 
needs of the applied arts to be fulfilled, but also 
there must be clarification as regards the flow pro- 
cesses during unsteady flows. 
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3. Cylinder in a Field of Sinusoidal Motion 

Forces on a cylinder admit an easier representation 
when the undisturbed portion of the flow, infinite 
in extent, is varying harmonically. Let the velocity 
be given by 

U= -U m cosat, (17) 

where U m is the semiamplitude of the current, T 
the period of the alternations, and <t=2t/T. The 
force on the cylinder per unit length F is in general 



F--=f(t,T,U m ,D, P ,v). 



(18) 



Grouping the variables on the basis of dimensional 
reasoning 

F Jt U n T U m D\ 

wr f \f' -w —y 



pUJ 
or introducing 

F 



o=2rt/T, 



pui 



J, U m T U m D\ 



(19) 
(20) 



where \J m D/v is a Reynolds number and U m T/D 
will be termed the "period parameter. " Bearing 
in mind that F is periodic, and that because of flow 
symmetry 

F(0)= -F(fi+T), 
we have 



pUID 



=A X sin 0+A 3 sin 30+ A 5 sin 50+. . 
+Bi cos 0+5 3 cos 30+B 5 cos 50+ 



(21) 



Here the coefficients A i} A 3 . . ., and B h B 3 . . . are 
independent of 0, and are at most functions of 
XJ m T/D and U m D/v. A sure method of approach 
in the analysis of the observed force curve is to re- 
sort to a Fourier analysis to determine the coeffi- 
cients Ai . 



and 



-B, 



2T F sin nd 
WW 

27r F cos nd 
pUID 



de 



de. 



(22) 



(23) 



Once the coefficients are obtained, their dependence 
on U m T/D and U m D/v may be established, provided 
the observational data are of sufficient number and of 
large extent. 

The above general and fundamental relation, eq 
(21), may be reconciled with eq (4), which is the form 
which Morison and coinvestigators Reid, Bret- 
schneider and others, have adopted in their numerous 
studies. Introducing U from eq (17) into eq (4) 



F 



pUW 4 



IT ~ 



Yj- sin — — |cos 0|cos 

U m L 



(24) 



By the rule of Fourier 

Jo 



cos 0|cos cos nddd 



|cos 0|cos fl=y^, 

71 = 



Jo 



where 



cos 2 nddd 
=a -\-ai cos 0+a 2 cos 20+a 3 cos 30+. 

a n =0 for n even, 



a n =(— 1)" 



n(n 2 — 4)t 



for n odd, 



di=— > a 3 =— z— j a 5 - 

07T l07T 



1057T 



(25) 



Introducing this in eq (21), and writing 

i?, 






B t 



(26) 






one has 

F 
TJ2T) =A l sin d+A 3 sin 30+ A B sin 50+ . , 

+B[\ cos 0| cos d+Bz cos 30+5; cos 50+ 



(27) 



Now eq (24) and (27) may be compared. One can 
write 



7r n Da_ A , A sin 30 , A sin 50 

i Vm • mr ^ z ^inT + ^ 5 ihT7 



and 



^=-Bl-Bl 



cos 30 



B' 5 cos 50 

cos 6\ cos I cos 0| cos 



or 



0,(0)=^ ¥f£ [A 1 +A Z +A 6 +2(A Z 



and 

C d (d)-- 



-2B[- 



+A 5 ) cos 20+2A 5 cos 40+ . . .] (28) 



[2(S;-B;)+4(5J-BJ)cos20 



cos 0| 



-45; cos 40+ . . .]. (29) 



Thus if A z , A 5 , and B' z , B' 5 vanish, the coefficients 
of mass and drag remain constant for all the phases 
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in the wave cycle and 



n _2_UmT A ..2 U M T C*" F sin OdO 
""^ D ^"tt 3 ' D Jo ptftO 



and 



a— «i— !X 



2ir F_cosJdd 
pUID ' 



(30) 



(31) 



In the event that these coefficients vary with the 
phase of the wave cycle, the values given by eq 
(30) and (31) are in a sense the weighted averages 



and 



G.-- (** C m (0)ain* MO (32) 

TTJO 

0*=+! f " ft(0) |cos 0| cos 2 A20. (33) 



With the above possibilities in mind, ii is prefer- 
able to adopt the expressions 

=A l sin0+#; cos |cos S\ + AR (34) 



or 



pU 2 m D 

F _* ' n Da . fl _Cd 
pU 2 m D~4 Cm 'U m smt) 2 ' 



Da C 

sin 0—^ |cos 0|cos d+AR, (34a) 



where A u B[, C mj and C d are constant, and AR has 
the value 

A/^=^4 3 sin 30+^s sin 50+^ cos 30+ /^ cos 50. (35) 

The function AR will be referred to as the remainder 
function, and then this remainder function is ob- 
tained by subtracting the computed values of A t 
sin0 and B[ |cos 0| cos from the observed F/pUmD. 
The remainder thus obtained may be examined in 
regard to its Fourier structure 4 and also as to its 
magnitude. 

4. Characteristics of the Experimental 
Waves 

The region under the nodal area of a standing wave 
that may be realized in a rectangular vessel furnishes 
a velocity field of simple harmonic motion in the 
velocity component U. This circumstance is not 
seriously modified even when the surges are moder- 
ately high. 

Taking the x-axis in the plane surface of the un- 
disturbed water, the z-axis vertical and upwards and 
the origin at one end of the basin, (see fig. 2), the 
surface elevation as reckoned from the undisturbed 
level, according to the second-approximation theory, 
from Miche [16], is 

ak 



h=a cos kx sino-tf+a — iVi cos 2kx~ 
a—r-N 2 cos 2kx cos 2crt, 



(36) 



where 
and 



N V 



N 2 = 



cosh 2kH 
= sinh 2kH 

cosh 2 kH (cosh 2kH+2) 



sinh 2 kH sinh kH 



Here k = T/L, L being the length of the basin; a= 
2x/T, J 7 being the period of oscillation; H the depth 
of water; and a the semiwave height, that is, the 
mean value of the extreme end deflections in a cycle. 
The expression for the period is the same as in the 
first-approximation theory, that is, 



a 2 =gk tanh kH. 



(37) 



Focusing attention on the basin end x=0, the surface 
displacement is 



h=asm<rt+a-T-Ni—a-^N 2 cos 2at; 



x=(). 



(38) 



Thus, the maximum elevation, occurring at£=7r/2<r, is 
h^a+a^lN.+m, (39) 

and the maximum depression, at £=3t/2<t, is 

The ratio of the elevation to the depression is 

|= -(i+^i^+m^jii-^iN.+m^y (40) 

and accordingly its value increases with wave height. 
The surface configuration for t=Q is 



h=(tj[N l —NtOOs2kx], 



t=0. 



(41) 
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Figure 2. Notation diagram for wave profile. 
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This represents a positive hump at the center of the 
basin and depressions at the ends. As a result, the 
duration of time that the surface of the water at one 
end of the basin is found to be above the undisturbed 
level is shorter than the duration that it is below. 
This matter has a bearing on the manner of fixing the 
reference time of the force cycles studied, and requires 
further discussion. 

At a small positive time t=T , the elevation h is 
nil, and this is the time when the wave in its upward 
surge reaches the undisturbed level. 

Since ar is a small angle, sin ar =aTo, and from 
eq (38) 

ar = (Nz-N^-kNz-m^Nz. (42) 



4 2 V 



16 



At a later time, t= 2 T /2 + r 1 , once more h = 0. This is 
the time when the wave in its downward surge reaches 
the undisturbed level. Since ar x is also a small angle, 
sin 0Ti=<rn, and from eq (38) 



and, thus, 



Tl = 



"TO 



(43) 



Let T t denote the duration of time that the surface 
of the water at the end of the basin, x=0, is above 
the undisturbed level, and T the duration below the 
same level. Accordingly, 



and 



By definition 



T t +T =T, 



2(T -T t ) m 



HP> 



T i = ^+T 1 ~T , 



and in view of eq (43) 



or 



and, thus 



T 



2T t _ 2ar 

T ~ TT 



2(T ~T i ) _ 4(7T 

T TT " 



Introducing the value of ar from eq (42) 

(44) 



If the instant when the upsurging wave at the end, 
x—0, reaches the level of the undisturbed water is 
observed, this then determines the instant t=T ; 
Since 4:t =T — T u the value of r may be obtained 
from the time durations that the water surface is 
below or above the still level. If on the other hand 
these observations have not been made, then r must 
be obtained from eq (42), introducing in it the wave 
height a of the observed surge deflections. 

The expressions for the particle velocities within 
the order of the approximations considered are from 
Miche [16], 



u= 



gak coshfc^+fiO 
a cosh kH 



sin kx cos at 



3ga 2 k 2 cosh2k(z+H) 
4 a sinh 2 kH sinh 2kH 



sin 2kx sin 2at (45) 



and 



oak sinh k(z+H) 7 
w=— — i 7 TT — cos kx cos at 

a cosh kH 

+7 - ^ 2 - S S h /^ ( -t?ir7 cos 2kx sin 2at. (46) 
4 a smh 2 kH smh 2k H 

At the vertical plane through the midsection of the 
basin, that is, at the plane x=L/2 or kx=ir/2, the 
velocities are 



oak cosM (z-\-H) 
u=-— i y tt cos at 

a cosh kH 



(47) 



and 



w=—--a 2 k 
4 a 



, sinh 2k(z+H) 
' sinh 2 kH sinh 2k H 



sin 2at. (48) 



Thus at the channel midsection, the horizontal 
component of the particle velocities is simple 
harmonic. The vertical component is also simple 
harmonic except that the frequency is twice as 
large. The effect of vertical velocity decreases with 
wave height. It is further reduced by lowering the 
object in the basin. Denoting the position of the 
object in the basin by Z\ and putting 



77 _ffflfc cosh k(zi+H) 
Um ~~a cosh kH ' 

the velocity components are 

u= — U m cos at 



(49) 



(50) 



and 



w=— - kH 



sinh k(Zi~\-H) a 



sinh 3 



ttt tt U m sin 2at. (51) 



It is inferred that w becomes less significant when 
kH is larger than 0.9. This limits the length of the 
basin for a given depth of water. For studies of 
wave forces in basins of greater length, or with water 
of less depth the present theory proves inadequate. 

All the experiments discussed subsequently were 
made in a basin of length i=242 cm and water 
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depth H=70 cm. The objects, cylinders or plates, 
were placed 25 cm below the water surface, that is 
^1= — 25 cm in the midsection plane ^=121 cm. 

For this condition, //&= 0.908, iVi = 1.054, and 
N 2 =3.322. From eq (40), the ratio of end deflections 
reduces to 



h 2 ~ 1—0.993 a/H 



(52) 



The graph of this equation is shown in figure 3, and 
values from observation are given by circles. The 
agreement between theory and observation is satis- 
factory for a/H less than 0.3. With this* restriction 
in mind, the value of the semiwave height a may be 
inferred from (39), that is, 



h 



= 1+0.993 a/#. 



(53) 



During the tests the elevation h\ was most easily 
observed. 

Fromeq (49) the relation between the current semi- 
amplitude and the wave height, in cgs units, is 



£7^3.43 a. 



(54) 



At z { the horizontal velocity is not uniform in the 
vertical direction. In the absence of a cylinder, 
with z measured in centimeters, 



1 dU m 
U m dz 



=0.00685. 



Thus, if A(J m be the difference in the maximum 
velocities at two points differing in elevation by D, 
then 

4^=0.48 D/H. 

For the largest cylinder used in the experiments, 
D=7.62 cm, the value of the ratio AUJU m is 0.052. 
The maximum value, during the cycle, of the 
vertical velocity component is given by 

w m /U m =0.3S a/H. 

The majority of the experiments were made with a 
less than 10 cm. For these cases, w m /U m is less 
than 0.055. 

From eq (44) the proportion of time that the sur- 
face of the water at one end of the channel is above 
or below the undisturbed level is given by 



^^ = 0.328—0.254 



(*)"■ 



(^) 



The graph of this expression is shown in figure 4 and 
values from observation are given by circles. For 
the observations, there was introduced into the 
basin at each end a parallel-wire resistance electrode, 
the bare parts of the wires being about 5 cm long 
and placed in a horizontal position just touching the 



surface of water at rest. The time that the electric 
current was traversing the electrodes gave the time 
that the water surface was above the undisturbed 
level, as in figure 5. 
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Figure 3. Variation of end deflections loith wave height. 
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Figure 4. Difference in the durations of the end elevations and 
depressions. 




Figure 5. Time record of the durations of the end deflections. 



429 



5. Dynamometer and Calibration 

The sketch of the dynamometer assembly is shown 
in figure 6. The rigid and massive base A for sup- 
porting purposes is firmly attached to the steel frame 
of the rocking basin directly above the water. The 
dynamometer itself consists of a pendular frame to 
which is attached the object to be immersed in water, 
a cylinder or a rectangular plate. The frame is 
constructed of brass angles and is strong enough to 
resist torsional and flexural deformation. The pivot 
depressions, located at the upper corners of the 
frame, consist of small bores of 1 mm in diameter in 
a bronze bedding. The bores are centered about 
polished steel conical points emerging from the 
supporting base. At a lower level two duraluminum 
annular rings of rectangular cross section are clamped 
to the frame and to the base. 

These rings constitute elastic elements for measur- 
ing the forces. To indicate the ring deformations 
two pairs of strain gages, SR-4, 120 ohms, are glued 
to each of the rings, inside and outside, and at dia- 
metrically opposite points. The four strain gage 
elements form the bridge which is led to a universal 
analyser. The latter is relayed to one of the chan- 
nels of a two-channel magnetic oscillograph. The 
second channel is reserved for timing observations. 
A similar connection is adopted for the other ring. 
By having four strain gages on each ring the sensi- 
tivity is increased and no corrections are needed for 
temperature changes. Two different sets of rings 
are used for measuring forces of different magni- 
tudes. The method of calibration may be inferred 
from the sketch in figure 7. The sum of the forces 
on the two rings equals 0.625 times the load applied 
to the frame. The ring deformations are examined 
for loads producing tension and compression. The 
indications of the ring deformations as read from the 
oscillograph record are linear as shown in figure 7. 
The calibrations were repeated before each run to 



guard against accidental changes in the strain-gage 
behavior. 
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70.0 



Figure 6. Dynamometer assembly (dimensions in centimeters). 




Figure 7. Calibration of the strain gages. 
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Figure 8. An example of oscillograph record of forces. 
Run 82, UmTJD =15,6. 
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6. Record of Forces and the Reduction 

An example of two oscillograph records of the 
forces, one from each ring, and of the timing is shown 
In figure S. The nearly sinusoidal traces relate to 
the forces acting on the rings; the others, in steps, 
give the time sequence. The incidence of the larger 
deflections indicates the time that the parallel wire 
electrode at the basin end x=0 was immersed; and 
the incidence of no deflection indicates the time that 
the electrode was out of the water. The point P 
where the greater deflections appear to commence 
gives the instant that the upsurging water reached 
the undisturbed level. Thus the point P gives the 
time £=t , the value of which was computed from 
eq (55), 4r = T — T i) after introducing the semiwave 
height of the wave. 'Phis value was transferred to 
the record to mark the origin of time, t=0, shown by 
the line AA\ The line BB' indicates the end of the 
wave cycle and corresponds to t= T. To establish 
the correspondence of the records from the two rings, 
the timing marks appearing at the lower edges of 
the records were used. 

At the time the record of the forces was being 
taken the wave elevation hi was read visually 
against a paper scale attached to the end wall of the 
basin. The water surface was readily discernable 
through the Incite walls of the basin. The magni- 
tude of the semiwave height a was deduced from 
eq (53), using the observed value //,. Maximum 
current V m was deduced from eq (54). 

The sum of the corresponding readings of the 
sinusoidal tracings in figure 8 gives the magnitude of 
the forces acting on the two rings when the calibration 
is applied. Taking moments about the dynamom- 
eter pivot point, the total force X on the cylinder is 
obtained. This is divided by the length of the cyl- 
inder to give F. The time history of the reduced 
force, F/pU m 2 D, is shown in figure 9. 
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Figure 9. .1// example of a curve of reduced forces on a cylinder. 
Run 82, U n T t Z>*15.6. 



7. Inertia and the Drag Coefficients of 
Cylinders and Plates 

Considering the force data in dimensionless form. 

such as shown in figure 9, the coefficients .1, and /V 
of eq (34) were determined by the method implied in 
eq (30) and (31). The desired integrations were 
carried out in the form of summations by giving to 
the differential multiplier dB the incremental value 
A0=O.O57r. The values of A x and J?/ thus found are 
entered in table 5 for the cylinders and in table 6 for 
the plates. Tables 1 and 2 contain the diameters of 
the cylinders, or the width of the plates, the values 
of the maximum currents and the water tempera- 
tures. Next the values of the inertia coefficient, 
C my were determined on the basis of eq (30), and the 
values of the drag coefficient, C dj on the basis of eq 
(31). These results are entered in table 3 for the 
cylinders and in table 4 for the plates. These tables 
also contain the Reynolds number U m D/i> and the 
period parameter U m T/I). 

The inertia coefficient C m varies from the theo- 
retical value 2 as the diameter of the cylinder is 
changed, or with a given cylinder as the maximum 
current is varied. Similar variations occur in the 
drag coefficient C[ h the changes being in the form 
of additions to the value experienced in steady flow. 
A correlation between the coefficients and Reynolds 
number U m D/v does not appear to exist. On the 
other hand, when these coefficients are related to the 
period parameter U m T/D, definite and regular de- 
pendencies arc 4 found. This is illustrated in figures 
10 and 11 for the cylinders, and in figures 12 and 13 
for the plates. 

Table 1. Cylinders 
[T-2.075secl 



Run 



1.. 
2__ 
3.. 
4__ 

5. 

6.. 
7_. 
8.. 
9.. 
lu- 
ll. 
12. 
13- 
14. 
15. 

16- 
17- 
18. 
19- 
20. 

21. 
22. 
23. 

24_ 

25. 

26. 
27- 
28- 
29- 



3 
3 
3 
3 
3 

3 
3 
3 
3 
2.5 

2.5 
2.5 
2.5 
2.5 
2 

2 
2 
2 
2 
1.5 

1.5 
1.5 
1.5 
1.5 

1.25 

1.25 
1.25 
1. 25 
1.25 



U n 



cm .vc 
36.2 
30.2 
27.7 
24.5 
21.1 

19.2 
15.8 
13.1 
10.0 
33.1 

27.4 
20.7 
13. 
10.3 
41.5 

35.4 
27.5 
19.1 
23.5 
53.2 

43.4 
33.4 

25.7 
19.4 
62.9 

54.5 
43.8 
35.7 
27.1 



C 
23.0 
22.0 
23. 
22.5 
22.5 

22.0 
22.0 
23.0 
23.0 
24.0 

24.0 
24.0 
24.0 
24.0 
24.0 

24.0 
24.0 

24.8 
24.8 
26.0 

26. 
26.0 
26.0 
26.0 
24.0 

28.0 
24.0 
29.2 
29.2 



Run 



30. 
31 
32 
33. 
34. 

35. 
36. 
37. 
38. 
39. 

40. 
41. 
78. 
79. 
80. 

81. 
82 
83. 

84. 
85. 

86. 
87. 
88 
89 
90 

91 
92 



in. 

l 
l 
l 

l 

0. 75 

.75 
.75 
.75 
.75 
.5 

.5 

.5 

1.75 
1.75 
1.75 

1.75 

1.5 

1.5 

1.5 

1.5 

0.5 
.5 
.5 
.5 
.75 

.75 
.75 
.75 



cm/ sec 
71.7 
58.7 
45.3 
36.0 
70.7 

63.8 
53.5 
45.3 
38.1 
73.4 

58.7 
48.0 
27.6 
24.0 
17.7 

14.4 
28.7 
25.2 
20.2 
14.6 

66.4 
54.8 
44.6 
32.6 
54.0 

49.6 
46.0 
41.0 



C 

30.0 
30. i) 
30.0 
30.0 
30.0 

30. 
30. 
30.0 
30. 
30. 

30. 
30.0 
22. C 
22.0 
22.0 

22.0 
20.5 
20.5 
20.5 

20.5 

21.0 
21.0 
21.0 
21.0 
12.0 

12.0 
12.0 
12. 



461441- 
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Table 2. Plates 

[T= 2.075 sec] 



Run 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


D 


U m 


e 


Run 


D 


Um 


e 


in. 

3 

3 

3 

3 

3 

2.5 
2.5 
2.5 
2.5 
2.5 

2 
2 
2 
2 
2 

1.5 
1.5 
1.5 


cm/sec 
14.1 
12.7 

9.8 
8.0 
6.4 

18.5 
15.8 
13.0 
10.3 
6.5 

21.6 
18.9 
16.1 
13.0 
9.9 

30.0 
25.0 
18.9 


° c 

23.0 
23.0 
23.0 
23.0 
23.0 

24.0 
24.0 
24.0 
24.0 
24.0 

27.0 
27.0 
27.0 
27.0 
27.0 

24.8 
24.8 
24.8 


60 

61 

62 

63 _ 


in. 

1.5 

1.25 

1.25 

1.25 

1.25 

1 
1 
1 

1 
0.75 

.75 

.75 

.75 

.5 

.5 

.5 
.5 
.5 


cm/sec 
13.7 
37.0 
29.5 
22.0 
16.1 

41.2 
34.5 
27.4 
18.9 
57.1 

47.0 
37.7 
27.4 
72.4 
63.6 

54.0 
45.3 

35.8 


°c 

24.8 
28.5 
28.5 
28.5 

28.5 

30.4 
30.4 
30.4 
30.4 
30.8 

30.8 
30.8 
30.8 
30.0 
30.0 

30.0 
30.0 
30.0 


64. 


65 

66 


67 


68 


69 


70 


71 


72 


73 


74 


75 


76. ._ 


77 





For the cylinders, as one passes from the small 
values of the period parameter to the larger values, 
the inertia coefficient commences to fall from the 
initial value 2 to a minimum value of 1.00 at 
U m T/D=15 and then graduallv increases to a value 
of 2.5 at U m TID=120. In regard to the coefficient 
of drag, there is an increase from the initial value 
0.9 to a maximum value 2.5 at U m TjD—\5 and then 
there is a gradual decrease to the value obtained in 
steady flow. It appears that for the cylinders the 
narrow region around U m T/D=l5 is a critical one. 

As regards the plates, the course of the variations 
of C m and C d with the period parameter is of a very 
different kind. It will be noticed that C m first in- 
creases, then decreases and finally rises again to a 
value of nearly 4.5. The most remarkable behavior, 
however, is in regard to C d . The coefficient of 
drag, starting with an unusually large value, 10, 
decreases rapidly at first and then gradually for 



Table 3. Cylinders 



Run 



1_. 
2.. 
3.. 

4_. 
5_. 

6.. 

7.. 
8.. 
9_. 
10 

11 
12. 
13. 
14. 
15 

16. 
17. 
18. 
19. 
20. 

21. 
22. 
23. 
24. 
25. 

26. 
27. 

28. 
29- 



C m 


c d 


U m T/D 


1.11 


1.24 


9.9 


1.44 


1.14 


8.2 


1.49 


1.32 


7.6 


1.70 


1.13 


6.7 


1.88 


1.00 


5.7 


1.95 


0.91 


5.2 


2.05 


1.23 


4.3 


2.10 


1.01 


3.6 


2.14 


0.70 


2.7 


0.74 


1.69 


10.8 


1.14 


1.61 


8.9 


1.71 


1.36 


6.8 


2.02 


1.15 


4.3 


2.06 


1.12 


3.4 


0.72 


1.73 


17.0 


.70 


1.98 


14.5 


.83 


2.18 


11.2 


1.50 


1.89 


7.8 


1.10 


1.97 


9.6 


1.02 


1.30 


29.0 


1.02 


1.49 


23.6 


0.85 


1.73 


18.2 


.74 


2.15 


14.0 


.87 


2.21 


10.6 


1.24 


1.15 


41.1 


1.27 


1.23 


35. 6 


1.40 


1.46 


28.6 


1.26 


1.52 


23.3 


0.87 


1.75 


17.7 



U m D/y 



293X102 

239 

225 

196 

169 

152 
125 
106 
81 
229 

189 
143 

90 

71 
230 

196 
152 
108 
132 
231 

188 
145 
112 
84 
218 

206 
152 
138 
105 



Run 



30 


1.66 


31.. _ 


1.70 


32 


J. 68 


33___ 


1.64 


43 


1.82 


35 ._ 


1.61 


36 .. 


1.63 


37 . 


1.64 


38 


1.84 


39 _ 


2.54 


40 


2.35 


41 . 


2.15 


78 


0.82 


79 


.84 


80 


1.41 


81— 


1.78 


82. ._ 


0.80 


83 .. 


.78 


84 


.87 


85 


1.46 


86. . 


2.52 


87. 


2.60 


88 . 


2.32 


89 .. 


2.26 


90 _. 


1.82 


91 


1.81 


92_._. 


1.81 


93 


1.76 



C d 



1.09 
1.29 
1.40 
1.49 
1.10 

1.19 
1.42 
1.45 
1.50 
1.07 

1.29 
1.42 
1.99 
2.08 
2.06 

1.75 
2.05 
2.28 
2.36 
2.18 

1.18 
1.31 
1.43 
1.54 

1.28 

1.38 
1.42 
1.54 



UmT/D 



58.6 
47.9 
37.0 
29.4 
77.0 

69.5 
58.3 
49.3 
41.5 
119.9 

95.8 
78.5 
12.9 
11.2 
8.3 

6.7 
15.6 
13.7 
11.0 

7.9 

108.4 
89.5 
72.9 
53.2 

58.8 

54.0 
50.1 

44.7 



U m D!v 



227X102 

185 

143 

114 

167 

151 
127 
107 
90 
116 

93 

76 
127 
111 

82 

67 
109 
96 

77 
56 

85 
71 
57 
42 
83 

76 
71 
63 











Table 4 


Plates 










Run 


C m 


C d 


U m T/D 


UmBlv 


Run 


Cm 


Cd 


U m T/D 


UmDjv 


42 


1.94 
1.74 
1.56 
1.51 
1.35 

2.28 
2.12 
2.00 
1.91 
1.57 

2.22 
2.44 
2.42 

2.17 
2.16 

0.95 
1.07 

2.08 


8.75 
8.81 
9.76 
10.21 
11.55 

5.50 
7.06 
8.01 
8.64 
11.44 

5.21 
5.48 
6.31 
7.25 
8.04 

4.11 

4.28 
4.61 


3.8 
3.5 
2.7 
2.2 
1.7 

6.1 
5.2 
4.2 
3.4 
2.1 

8.8 
7.7 
6.6 
5.3 
4.1 

16.3 
13.6 
10.3 


114X102 
103 

80 

65 

52 

128 
109 

90 

71 

45 

128 
112 
95 

77 
59 

127 
106 
80 


60 


2.51 
2.14 
1.07 
1.43 
2.25 

2.45 
2.10 
2.01 
1.56 
3.17 

2.88 
2.89 
2.71 
4.96 
4.09 

4.00 
3.58 
3.70 


5.15 
3.25 
3.94 
4.09 
4.43 

3.13 
3.55 
3.68 
4.38 
2.43 

2.86 
3.06 
3.36 
1.81 
2.03 

2.32 
2.45 
2.59 


7.5 
24.2 
19.3 
14.3 
10.5 

33. 6 
28.2 
22.4 
15.4 
62.2 

51.2 
41.1 
29.9 
118.2 
104.0 

88.3 
74.0 

58.6 


58X10 2 
142 
113 

84 

62 

131 
110 

87 

60 
138 

113 

91 

66 
114 
101 

85 

72 
57 


43 


61 

62 


44 


45 


63. _. 


46 


64 


47 


65... 


48 

49 


66._. 


67 

68 

69 


50 

51 

52 

53 

54 

55 

56 

57 


70 

71 

72 

73 

74 


75 

76 

77 


58 


59 
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Table 


5. Cylinders 






Run 


Ax 


B[ 


A* 


Bi 


sU 


B' b 


UmT/D 


x 


0.56 
.87 
.97 
1.26 
1.62 

1.84 
2.35 
2.91 
3.87 
0.34 

.63 
1.25 
2.34 
3.03 
0.21 

.24 
.36 
.95 
.57 
.17 

.21 
.23 
. 26 
.41 
.15 

.18 
.24 
.27 
.24 
.14 

.17 
.22 
.27 
.12 
.11 

.14 
.16 
.22 
.10 
.12 

.14 
.31 

.37 

-1 

L.30 

0.2B 

.28 
38 

.\)\ 
.11 

.14 
.111 
.21 
.15 
.17 

.18 
.19 


-0.62 
-.57 
-.66 
-.56 
-.51 

-.4(> 
-.62 
-.51 
-.35 

-.84 

-.SO 
-.(is 
-.58 
-.50 
-.96 

-.99 
-1.09 
-0.95 
-.99 
-.65 

-.74 

-.87 

-1.08 

-1.10 
-0.58 

-.62 
-.73 
-.76 
-.87 
-.55 

-.64 
-.70 
-.74 
-.55 
-.60 

-.71 
-.72 
-.75 
-.54 
-.65 

-.71 
-.99 
- 1 . 04 
-1.03 
-0.88 

-1.02 
-1.14 
-1.18 
-1.09 
-0. 59 

-.05 
-.72 
-.77 
-.1)4 
-.69 

-.71 
-.77 


-0. 11 
-.00 
-.OS 
-.03 
-.02 

+.02 
.01 
AY.\ 
.01 

-.19 

-.14 
-.07 
+ .10 
.01 
-.18 

-.22 
-. 19 
-.14 
-.15 
-.05 

-.07 
-. 15 
-.22 
-.20 
.00 

-.01 
-.01 
-.04 
-.13 
+.02 

.02 
.01 
.01 
.03 
.01 

.01 
.(11 
.02 
.03 
.02 

.01 
-.18 
-.19 
-.15 
-.07 

-.22 
-.24 
-.23 
-.15 
+ .03 

.(12 
.01 
.01 
.02 
. 02 

.02 
.02 


-0.04 
-.03 

.00 

+ .05 

.05 

.06 

.12 
. 14 

.08 
-.02 

-.05 

-.02 

+.13 

.14 

.00 

-.01 
-.05 
-.04 
-.06 
.00 

.00 
-.05 
-.04 
-.07 
+ .01 

.01 

. 02 

.00 

-.04 

+ .01 

.01 
. 03 
.02 
.01 
.01 

.02 

.02 
.02 
.00 
.01 

.02 
-.05 
-.02 
-. 10 
-.10 

-.03 

-.01 
-.05 

Ml 

+ .01 

.01 
. 02 
.04 
.01 
.01 

. 02 
.02 


-0.03 
-.02 
-.02 

-.01 
.00 

.00 
-.01 
+ .01 
-.02 
-.04 

-.04 
.00 

+.01 
.00 

-.04 

-.02 
-.02 
-.02 
-.01 
-.01 

.(H) 
-.04 
.00 
.00 
.00 

.(H) 
.(H) 
.00 
-.04 
.00 

-.01 

.00 

.00 

-.01 

-.01 

.00 
-.01 
.00 
.00 
.00 

.00 
-.01 
-.03 

-.01 
.00 

-.04 

-.01 

-.02 

.00 

.110 

-.01 
-.01 

+ .01 
.(M) 
.00 

.00 
.00 


-0.08 
-.06 
-.04 
-.03 
-.03 

-.01 

-.02 

.00 

+ .02 
-.08 

-.09 
-.06 
-.01 
+.01 
-.05 

-.07 
-.12 
-.11 
-.11 
-.03 

-.04 
-.07 
-.08 
-.11 
-.01 

-.02 
-.02 
-.05 
-.07 
-.01 

-.01 
-.02 
-.02 
-.01 
-.01 

-.02 
-.01 
-.01 
-.01 
-.01 

-.01 
-.10 
-.12 
-.13 
-.10 

-.06 
-.10 
-.14 
-.12 

-.01 

-.01 
-.01 
-.01 
-.01 
-.01 

.00 
-.01 


9.9 
8.2 
7.6 
6.7 
5.7 

5.2 
4.3 
3.6 
2.7 
10.8 

8.9 
6.8 
4.3 
3.4 
17.0 

14.5 
11.2 

7.8 
9.6 
29. 

23. 
18.2 

11.0 
10. 
41.1 

:-jf>. (i 

2S. (i 
23.3 
17.7 
58. 6 

47.9 
37.0 
29.4 
77.0 
69.5 

58.3 
49.3 

M.:. 
119.9 

Of), s 

78.5 
12.9 
11.2 
8.3 
6.7 

15. 
13.7 
11.0 
7.9 
108.4 

89. 5 
72. 9 
53. 2 
58. 5 
54.0 

50. 1 
44.7 


2 

3 

4 


5 

6 


7 

8— 


9... 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 . 


27 

28- 


29 

30- . 


31.. - 


32 

33 

34 

35 . 


36- 


37 . 


38 . 


39.. 


40- 


41 


78- 


79_. 


80- 


81- 


82- . 


83 


84 


85 


86— 


87- _ 


88.. 


89- 


90 

91 — 


92.. 


93.. 





\?~ 



Table 6. Plates 



Run 



42 

43 

44 

45 

46— 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66— 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 



2.50 
2.49 
2.86 
3.42 
3.83 

1.86 
2.03 
2.32 
2.81 
3.63 

1.24 
1.56 
I s! 
2.01 
2.63 

0.29 
.39 
1.00 
1.66 
0.44 

.27 
.49 
1.05 
0. 30 
.37 

.44 
.50 
.25 
.28 
35 

.45 
.21 
.19 
.22 
.24 



B[ 



-4.38 
-4.41 
-4.88 
-5.11 

-5.77 

-2.75 
-3.53 
-4.01 
-4.32 
-5.72 

-2.61 
-2.74 
-3. 16 
-3.63 
-4.02 

-2.05 
-2.14 
-2.31 

-2. 58 

- 1 . 62 

-1.97 
-2. 04 
-2. 22 
-1.50 
-1.77 

-1.84 
-2.19 
-1.21 

- 1 . 43 
-1.53 

-1.08 
-0. 91 
-1.01 
-1. 10 
-1.22 



Ai 



II :-,S 

.58 
.47 
.53 
.43 

.37 

.45 
.53 
.68 
.61 

.35 
.31 

.38 
.45 
. 66 

.29 
.30 

.29 

.no 
. n 

is 
.HI 
.HO 
.08 
.00 

. 10 
.32 
.05 
.00 

.os 

. II 
.00 

.06 

. 05 
.05 



0.56 
.47 
.40 
.45 



.58 
.59 
.64 
.52 
.55 

.42 

.49 
.56 

.54 



.33 
.31 
.43 

.54 
.18 

. 19 
. HO 
.45 
. 13 
. 10 

. 1!) 
.38 
.05 

(is 

.11 

.13 
.03 
.03 

. 05 
.05 



A& 



o.oi 
.04 
-.02 
-.01 
-.02 

-.08 
-.04 
-.02 

+.07 
.01 

-.05 
-.00 
-.05 
-.01 
.00 

+.02 
.00 
-.08 
-.07 
-.05 

+ .02 
-.02 
-.00 
-.05 
-.04 

-.02 
-.01 
-.02 
-.OH 
-.04 

-.03 
-.01 
-.01 
-.02 
-.02 

-.03 



B' 5 



-0. 05 
-.11 
-.18 
-. 17 
-.10 

-.03 
-.02 
-.02 
-.00 
-.14 

+.02 
-.03 
-.02 
-.03 
-.01 

+.07 

.03 

. 00 

-.04 

-.02 

+.07 

.01 

-.01 
-.01 
-.01 

+.03 
06 
-.04 
-.04 
-.04 

-.01 
-.03 
-.03 
-.03 
-.02 

-.03 



U m T/D 



3.8 
3.5 
2.7 
2.2 

1.7 

6.1 
5.2 
4.2 
3.4 

2.1 

8.8 
7.7 
0. (i 
5.3 
4.1 

16.3 

13.6 

10. H 

7.5 

24.2 

[9.3 
1 1.8 
10.5 
33.6 

28.2 

22. I 
15.4 
62. 2 
51.2 
41.1 

29. 9 

US. 2 
101.0 

ss. 3 
74.0 

58.6 
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Figure 10. Variation of inertia coefficient of cylinders. 
Diameter (inches) : 3 2.5 2 1.75 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D O • A ■ ♦ + 



U m T/ D 

Figure 11. Variation of drag coefficient of cylinders. 
Diameter (inches): 3 2.5 2 1.75 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D O • A ■ ♦ + 

increasing period parameter. The final value is 
almost identical with that found for steady flow. It 
is perhaps important to mention thai O'Brien 
and Morison [17] noted equally large values of drag 
coefficients for spheres subjected to t he action of pro- 
gressive waves. It will be noted that the larger 
values of C d are associated with the smaller values of 
C m and the larger values of C m with the smaller 
values of C d . Because the drag coefficient is large 
when U m TjD is small and the variation of C m is 
relatively moderate, the wave forces on plates are 
essentially due to drag, and the inertia effects play 
a small role almost independent of the period 
parameter. 
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Figure 12. Variation of inertia coefficient of plate 

Diameter (inches): 3 2.5 2 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D # ▲ ■ ♦ + 






U m T/D 



Figure 13. Variation of drag coefficient of plates. 

Diameter (inches): 3 2.5 2 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D • A ■ ♦ + 



8. Variations of the Remainder Function 
and C m and C d During the Wave Cycle 

The values of C m and C d given in tables 3 and 4 are 
average values for the entire wave cycle, and in some 
cases local values may differ from the average. 
Where the inertia and drag coefficients, C m and C d , 
each have the same constant value at all phases of the 
wave cycle, eq (24) should suffice to describe ade- 
quately the magnitude of the forces at every phase. 
On the other hand, should C m and C d vary with the 
different phases, the forces are better represented by 
eq (34a). The variations in C m and C d should lead 
to the remainder force function AR. The examples 
of the remainder function AR are given in figures 
14 and 15 where AR is the difference between 
f IpUmD—Ai sin and B[ cos 6 | cos 6 | in conformity 
with eq (34). Once a curve of AR as a function of 6 
is obtained, its structure in Fourier components may 
be considered and the coefficients A 3 , A 5 . . ., and 
B'z, B' b may be obtained. These determinations 
are given in tables 5 and 6 and in figures 16, 17, 18, 
and 19. 

Now for the determinations of the local values of 
C m and C d , two methods are available. The first 
gives the point values of the coefficients in a wave 
cycle as determined from the observed values of 
F/pUmD, using eq (24). Two sets of evaluations 
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Figure 14. Evaluation of remainder force AR for a cylinder. 

Run 82, U m T/D = 15.Q. 
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Figure 15. Evaluation of remainder force AR for a plate. 
Run 54, (7 m T/i>=6.6. 
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Figure 16. Variation of coefficients of the remainder force of 
cylinders. 

Diameter (inches): 3 2.5 2 1.75 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D O • Alt + 
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Figure 17. Variation of coefficients of the remainder force cf 
plates. 

Diameter (inches): 3 2.5 2 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D # ▲ ■ ♦ + 
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Figure 18. Variation of coefficients of the remainder force of 
cylinders. 

Diameter (inches): 3 2.5 2 1.75 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D O # A ■ ♦ + 
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Figure 19. Variation of coefficients of the remainder force cf 
plates. 

Diameter (inches): 3 2.5 2 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D • A ■ ♦ + 

were made, the basic suppositions being as follows: 
It was assumed in the first evaluation that for 
d 1 = Tr/2 J ra and 2 =W2 — «, whore a is an angle less 
than 7r/2, the coefficients C m and (\, each have equal 
values, since these are the phases where the accelera- 
tions, du/dt, arc equal, and the currents u, are equal 
in absolute value although of opposing signs. This 
is true also for 6i = 3ir/2 + a and 2 =33r/2 — a. In 
the second evaluation, it was assumed that for 
= tt + P and d=w—l3, where £ is an angle less than w, 



the coefficients C m and C d each have equal values, 
since these are the phases where the currents, >/, are 
equal and the accelerations, du/dt, are equal in 
absolute value, although of opposing signs. Also 
since we know the values of the coefficients Ai, A Zl 
A 5 and B[, B' z , B' 5 , the curves of C m and (\, as 
function of may be obtained by using eq (28) and 
(29). The latter is the second method and is mathe- 
matically equivalent to the assumptions made above. 

In the cylinder data the agreement between the 
observed and computed force is satisfactory when 
U m TjD is small. The computation was based on oq 
(24), introducing the values of C m and C d from the 
tables. Figure 20 illustrates this agreement. The 
local values of C m (d) and C d (0) for this case are shown 
in figure 21. The first determinations discussed 
above are shown by circles and squares, whereas the 
curves are determined by the second method. It is 
seen that C m (d) is independent of the phase t/T and 
that the coefficient C d (6) is constant except in short 
ranges of the phases t/7 1 =0.25 and 0.75. This is 
expected, for at these phases the current u vanishes. 
The values of C m and C d determined by eq (30) and 
(31) are given in the caption. 

The agreement between the observed and com- 
puted forces is also satisfactory when the period 
parameter is large. This is illustrated in figure 22. 
The local values of the coefficients for this ease are 
shown in figure 2\\. Here again, allowing small devi- 
ations, (\ n (6) is practically independent of the phase 
t/T and differs very little from the value given in 
table 3. On the other hand, considerable variations 
are obtained between the observed and. computed 
values of the forces in those cases where the period 
parameter is near U m T/D=\5, as shown in figure 24. 
The local values of the coefficients for this case are 
shown in figure 25. Now C m (6) varies considerably 
with the phase t/T, the smaller values occurring at 
t/T=0.0, 0.50, and 1, and the larger values at 
t/T=0.25 and 0.75. Also, C d {6) appears to be con- 
siderably augmented at the phases where the velocity 
u vanishes, that is, at t/T=Q.25 and 0.75. The ex- 
ample shown is typical for all the cases where 
U m T/D is in the neighborhood of U m T/D=15. In 
the example shown in figure 25, C m (6) shows large 
negative values at the points tjT—0, 0.5, and 1.0. 
The significance of this is not clear. It is believed, 
however, that the presence of negative values is not 
related altogether to the observational methods that 
were used. 

For the plates deviations were always found be- 
tween the observed values of the forces and the values 
computed on the basis of eq (24). An example is 
given in figure 26. The local values of C m (B) and 
C d (6) for this case are shown in figure 27. 
An additional example is given in figure 28. What 
is shown in these figures is typical for all the runs 
made with the plates. The coefficient (\ n {6) under- 
goes considerable variation in value for varying t/T, 
the greater values occurring at t/T=0, 0.5, and 1.0 
and the lesser values at t/T=0.25 and 0.75. Further- 
more, the increase in C d (B) at the points tlT=0.25 
and 0.75 is very decided. 
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Figure 20. Comparison of measured and computed forces on a 
cylinder. 
Run 9, U m T/D= 3.0. 
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Figure 21. ^.w example of variation of the inertia and drag 
coefficients of a cylinder during a wave cycle. 
Run 9, UmT/D=3.0, C m =2.U, C<*=0.70. 
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Figure 22. Comparison of measured and computed forces on a 
cylinder. 

Run 93, UmT/D=4A.7. 
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Figure 23. An example of variation of the inertia and drag 
coefficients of a cylinder during a wave cycle. 
Run 93, UmT/D=U.7, C»=1.76, C d =1.54. 




t/T 



Figure 24. Comparison cf measured and computed forces on a 
cylinder. 
Run 82, VmTJD =15.6. 
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Figure 25. An example of variation of the inertia and drag 
coefficients of a cylinder during a wave cycle. 

Run 82, U m T/D = 15.6, C m =0.80, C d =2.05. 

For the cylinder data, as long as the period parame- 
ter is sufficiently small, or sufficiently large, the forces 
may be computed on the basis of eq (24); the re- 
mainder function, AR, is small. For period parame- 
ters in the neighborhood of the critical value, 
U m T/D=15, the representation of forces is more 
exact using eq (34a); the remainder function is of 
significance. For the plate data the remainder may 
not be disregarded, in particular when the period 
parameter is small. 
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Figure 26. Comparison of measured and computed forces on a 

plate. 

Rim 54, U m T/D=QA}. 
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Figure 27. An example of variation of the inertia and drag 
coefficients of a plate during a wave cycle. 

Run 54, U m T/D = 6.G, C TO =2.42, C d =6.31. 



c m (0) 



CA9) 










o 


D O 

o 


o 

o 


o 








o o 


\ 


^-9- 






D 




^8-6' 


5/0 



□ 























-9— e— e— e— ©— ©- 



J L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

t /T 

Figure 28. An example of variation of the inertia and drag 
coefficients of a plate during a wave cycle. 

Run 69, U m T/D=62.2, C m =S.17, Cd=2.43. 



9. Flow Pattern Around Cylinders and Plates 

The flow patterns around the cylinders and plates 
for varying values of U m T/D were examined, because 
they may have had a bearing on the fact that the 

nature of forces during a cycle is significantly affected 
by the period parameter. The flow pattern was 
visually examined by introducing a jet of colored 
liquid 011 one side of the immersed object. The 
disposition of the liquid close to the object during 
the cyclic motion was recorded by a still camera and 
also by a motion-picture camera. Some of these 
pictures are shown in figures 29 and 30. 

Figure 29, a and b, were taken with the 3-inch 
cylinder, the first corresponding to U m T/D=4, the 
second to a larger value U m T/D=10. When the 
period parameter is small there is no separation of 
flow; the liquid near the cylinder clings to the cylin- 
der, and the partitioning of flow from above and 
below is symmetrical. It. will be remembered thai at 
low period parameter the inertia coefficient is about 
equal to the theoretical value 2, and drag is negligible. 
As U m T/D is increased there is separation of flow at 
the top surface of the cylinder during the relatively 
longer time that the flow continues in one direction, 
Although not. visible in the picture, somewhat later, 
but prior to the reversal of current, liquid coining 
around the cylinder from below moves upward and. 
although transforming into an eddy, remains close 
to the cylinder. 

Figure 29, c, illustrates the flow pattern for 
U m T/D=\7 with the 2-inch cylinder. Note the 
complete separation at the upper surface of the 
cylinder with the following flow around the lower 
surface directed upward with the subsequent eddy 
formation. 

A completely different picture is obtained for large 
period parameter, as shown in figure 4 29, d, taken 
with the K-inch cylinder, U m T/D=llO t Here one 
is confronted with the regular Karman vortices. 
The eddies are separated alternately from above 
and below. 

With plates the flow patterns are decidedly 
different, especially for small period parameter. 
Figure 30, a and b, show the 3-inch plate, the first 
corresponding to U m T/D=], the second to a larger 
value , U m T/D = 4 . Eddies a re formed almos t simul- 
taneously at the upper and lower edges of the plates. 
For the smaller value of U m T/D the eddies are 
apparently concentrated nearer the edges of the 
plate. Perhaps the large values of the drag coeffi- 
cient for small period parameter are associated with 
the behavior of the eddies in this case, but the ques- 
tion is left open for another occasion. 

Figure 30, c, illustrates the flow pattern for 
U m T/D=15 with the l#-inch plate, Here the eddy 
formation is no longer symmetrical, the separation 
occurring first at the upper edge of the plate followed 
by an eddy formed at the lower edge, remaining close 
to the plate. 

Again the Karman vortices are obtained for large 
period parameter as shown in figure 30d taken with 
the K-inch plate, U m T/D=110. 
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Figure 29. Flow patterns around cylinders. 

(A) Z>=3in., V m T\B=\\ (B) Z>=3 in., U m TID = 10; 
(C) D = 2in., U m T/D = 17; (D) Z>=0.5 in., [7 W T/Z> = 110. 
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cure 30, Flow patterns around plates. 

(A) D = 3in., U m T/D = l; (B) Z>=3 in., U m T/D = 4; 

(C) Z>=h5in., U m TjD=\b- (D) 2>=0.5 in., L T m T/Z) = 110. 
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The eddy appearances discussed above suggest the 
following interpretation as to the physical meaning 
of U m T/D. If one defines a length, /, as the distance 
that a fluid particle would move in one direction in 
the absence of the cylinder, l=U m T/ir. Thus, 



and accordingly the period parameter is proportional 
to the ratio of the distance traversed by a particle 
during a half cycle to the diameter of the cylinder. 
When the period parameter equals 15, l/D is 4.8. 
Perhaps when U m T/Di8 smaller than 15, the distance 
traveled by a particle is not large enough to form 
complete eddies. When it equals 15, the distance 
suffices to form a single vddy, and when much larger 
than 15 the greater distances allow the formation of 
numerous vortices of the Karman vortex street. 
One can hardly refrain from pointing to the similarity 
between (he period parameter and the Strouhal 
number, and as suggested by McNown and Keulegan 
[18], the product of Strouhal and period parameter 
numbers furnishes an alternate parameter as service- 
able as the period parameter number. If T s be the 
duration for the shedding of a single eddy, then the 
Strouhal number fD/U=S may be written as 
D/(2T S U) = S, since the number of alternative eddies 
shed during a second is 2/ and 2fT s equals 1 second. 
One may suppose that the relation is satisfied 
approximately also for sinusoidal motions, provided 
U is replaced by U m /2. Hence, the Strouhal number 
for sinusoidal motion is DI(U m T s )=S. 

Multiply the two sides by the wave parameter 
number U m T/D, 

T U T 

rr d m 

For cylinders, ignoring the dependence of S on the 
Reynolds number, 

U T 

As noted previously for the cylinders, C m attains its 
least value, slightly less than unity, at about 
U m T/D=12.5. This corresponds to the condition 
that T/T s = 2, nearly, and suggests that during a half 
cycle, that is, during a complete motion of fluid 
particle in one direction, a single eddy is formed and 
is separated (see also the figure 29, b). Obviously, 
the process of eddy shedding has a very significant 
bearing on the variations of the so-called coefficients 
of mass and drag, and account needs to be taken of 
this in the theoretical formulation of the basic 
process. 

10. Maximum Force During a Wave Cycle 

In engineering applications the main interest is in 
(he magnitude of the maximum force experienced 
during a wave cycle. If the remainder function is 
neglected, (he expression 



F/ P UID=A 1 sin d+B[\cos 0|cos 6 



instead of the eq (24), may be utilized to evaluate 
the maximum force F m and also i(s phase. If the 
maximum force FjpUffi occurs at 6=d m , the phase 
may be defined as 

The maximum value of the computed force is given 
by 

F 

pIPD = Al Sm ^ + ^ll C0S 0J COS Om, 



where 6 m satisfies the relation 



A l +2B[m\d m 



=O,orsin0 m = — ^±, 



for !<*„<!• 



As the coefficients A x and B'[ are functions of U m TjD 
only, then FJpUiD and <£ both are functions of 
U m T/D. For greater accuracy, the remainder 
function AR must be considered, but then the 
evaluations become somewhat involved. If these 
evaluations are made, the maximum force and phase 
are again functions of the period parameter. 

An alternative procedure is the direct establish- 
ment of the maximum force and phase by merelj 
taking these quantities from the reduced force 
curves of this investigation. Such readings for the 
cylinders are given in figure 31 and for the plates in 
figure 32. 
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Figure 31. Variations of the magnitude and phase of the 
maximum force on cylinders. 

Diameter (inches): 3 2.5 2 1.75 1.5 1.25 1 0.75 0.5 

Corresponding symbol: O A D <> • A ■ ♦ + 
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Figure 32. Variations of the magnitude and phase of the 
maximum force on plates. 

Diameter (inches) : 3 2. 5 2 1. 5 1. 25 1 0. 75 0. 5 

Corresponding symbol: O A D • A ■ ♦ 4- 

Table 7. Cylinders 





U m TjD 


F m 
pU 2 m D 


(degrees) 


U m TjD 


F m 


(degrees) 


pUiD 


2.5 


4.00 


82.0 


12.5 


1.28 


-6.8 


3.0 


3.45 


81.2 


15.0 


1.15 


-6.0 


3.5 


3.02 


80.8 


17.5 


1.03 


-5.0 


4.0 


2.65 


80.0 


20.0 


0.93 


-3.0 


4.5 


2.34 


79.0 


25.0 


.80 


1.0 


5.0 


2.10 


78.0 


30.0 


.73 


4.0 


5.5 


1.83 


76.6 


35.0 


.70 


6.5 


6.0 


1.60 


75.0 


40.0 


.68 


8.0 


6.5 


1.42 


72.5 


50.0 


.66 


8.0 


7.0 


1.30 


65.0 


60.0 


.65 


8.0 


7.5 


1.20 


35.0 


70.0 


.63 


8.3 


8.0 


1.20 


5.0 


80.0 


.63 


8.7 


9.0 


1.25 


-3.0 


90.0 


.62 


8.9 


10.0 


1.28 


-6.0 


100.0 


.62 


9.0 



Table 8. Plates 



UmTID 



2.0 
2.5 
3.0 
3.5 
4.0 

4.5 
5.0 
5.5 
6.0 
7.0 

8.0 
9.0 



F m 



6.50 
6.00 
5.55 
5.20 

4.85 

4.55 
4.30 
4.05 
3.82 
3.43 

3.10 
2.86 



(degrees) 



24.0 
25.6 
27.0 
28.2 
29.5 

30.8 
32.0 
33.0 
33.6 
34.9 

35.5 
35.5 



U m T/D 



10.0 
12.5 
15.0 
17.5 
20.0 

25.0 
30.0 
35.0 
40.0 
50.0 

60.0 
70.0 
80.0 
90.0 
100.0 



F m 
pU m D 



2.71 
2.44 
2.25 
2.10 

1.96 

1.76 
1.63 
1.53 
1.45 
1.33 

1.25 
1.18 
1.11 
1.06 
1.02 



(degrees) 



34.9 
32.0 
28.5 
26.0 
23.5 

18.3 
13.6 
11.1 
9.7 



9.0 
9.5 
10.1 
10.8 



For reference purposes, the data of the curves is 
given in tables 7 and 8, and can be used directly. 
In a future communication the forces on cylinders 
held in vertical positions will be computed on the 
basis of the data in these tables and will be com- 
pared with the laboratory observations already 
completed as a matter of concrete illustration. 



Washington, December 2, 1957. 
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